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ABSTRACT 



A method of defining a toleranced process based on a 
nominal process, applicable to any manufacturing process 
wherein the output is dependent on a process having a 
number of input factors which are subject to variation, and 
have a mean and standard deviation. The method comprises 
the steps of: representing the variability of the response of a 
system to the actual distribution of at least one of said factors 
and parameterized by at least one parameter; performing at 
least one step in the manufacturing process for a given lot of 
wafers; evaluating actual variability and nominal variability 
of said components against a schedule of parameters in an 
array; modeling the output of the evaluation to determine the 
manner in which to continue processing of the lot. 

32 Claims, 6 Drawing Sheets 
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LOT SPECIFIC PROCESS DESIGN (15) William D. Heavlin and Luigi Capodieci, "Calibra- 

METHODOLOGY tion and Computer Experiments," 1997 American Sta- 
tistical Association Proceedings, Section on Physical 
and Engineering Sciences, Anaheim, August 1997 

BACKGROUND OF THE INVENTION 5 (scheduled for publication Summer 1998), pp. 58-63. 

-* r** u c *i_ i (16) Heavlin, W. D., "Statistically Based Process 

1. Field of the Invention Windows, 1998 American Statistical Association 
The invention is directed to the mass production of Proceedings, Section on Physical and Engineering 

semiconductor devices. The invention is more particularly Sciences, pp. 216-221. 

directed to the problem of processing lots of wafers sharing 1Q K^^pic, Z., Heavlin, W. D., and Kyser, D., U.S. 

one or more common defects and to resolving those defects. Pat No 5)646f8 7o issued Jul. 8, 1997, entitled 

2. Cross Reference to Related Publications "Method for Setting and Adjusting Process Parameters 
The following publications are cited here for purposes of to Maintain Acceptable Critical Dimensions Across 

reference: Each Die of Mass-Produced Semiconductor Wafers." 

(1) S. Kaplan and L. Karklin, "Calibration of Lithography 15 (18) Krivokapic, Z., Heavlin, W. D., and Kyser, D., U.S. 
Simulator by Using Substitute Patterns," Proceedings Pat, No. 5,655,110 issued Aug. 5, 1997, entitled 
on Optical/Laser Microlithography VI, SPIE 1927, pp. "Method For Setting And Adjusting Process Param- 
847-858, 1993. eters to Maintain Acceptable Critical Dimensions 

(2) C. Mack and E. Charrier, "Yield Modeling for Across Each Die of Mass-produced Semiconductor 
Photolithography," Proceedings of OCG Microlithog- 20 Wafers." 

raphy Seminar, pp. 171-182, 1994. ^ ac ^ °f tne aforementioned publications is hereby incor- 

(3) TMA DEPICT, Two-Dimensional Process Simulation b y ref e rc nce 
Program for Deposition, Etching, and 3^ Description of the Related Art 
Photolithography, version 3.0, Technology Modeling Modern, high-density integrated circuit (IC) devices are 
Associates, Inc., Palo Alto, Calif., 1993. 25 mass-produced with large numbers of critically- 

„ iiT^r, . . - . . . dimensioned features. In manufacturing, it is desirable to 

(4) Mandel, J, The Statistical Analysis of Experimental ^ respectivc CTitical dimensions of each die 

Data, Chapter 12, Wiley, N.Y., 19o4. within a plurality of mass-produced IC wafers constrained to 

(5) Z. Krivokapic and W. D. Heavlin, "Predicting Manu- ceftain respeclive values in order to assure desired operating 
factunng Variabilities for Deep Micron Technologies: 3Q speeds and 0 p era tional characteristics of the produced IC. 
Integration of Process, Device, and Statistical Each f eature on eacD \q die in a mass-produced wafer is 
Simulations," in Simulation of Semiconductor Devices the product of a succe ssion of many process steps. Each 
and Processes, 5, S Selberherr, H Stippel and E process step is controlled by a combination of variable 
Strasser, eds, pp. 229-232, Springer- Verlag, N.Y., process parameters. 

1993. 35 Different combinations of variations in process parameter 

(6) W. D. Heavlin and G. P. Finnegan, "Dual Space can occur 0 n a random basis across the numerous process 
Algorithms for Designing Space-filling Experiments," stC ps of a mass-production line, on a die-by-die basis. This 
Interface 1994, Research Triangle, North Carolina, introduces noise into the uniformity of the product outflow 
June 1994. of the production line. Sometimes a specific combination of 

(7) B. D. Ripley, Spatial Statistics, pp. 44-75, Wiley, 40 process parameter deviations is relatively innocuous, some- 
N.Y., 1981. times it is not. It all depends on which process parameters 

(8) A. B.Owen, "Controlling Correlations in Latin Hyper- are deviated for a given IC die and how their respective 
cube Samples," Journal of the American Statistical process steps interrelate to establish critical dimensions on 
Association, Vol. 89, No. 428, pp.1517-1522, Decern- that given die. 

ber 1994. 45 One example of such multiple, interrelated process steps 

(9) W. D. Heavlin, "Variance Components and Computer thal ma y be useful ™ understanding the interaction of factors 
Experiments," 1994 ASA Proceedings, Section on in tne P resent application, are those typically employed to 
Physical and Engineering Sciences, Toronto, August define a P allera of conductive lines deposited across an 
1994 insulator of an IC chip. 

(10) A. R. Neureuther and F. H. Dill, "Photoresist Mod- 50 Firs '- a dieIe , c ( , lric la / er c of ' ' 8« era «y non-plinar form is 
eling and Device Fabrication Applications," Optical created ^ wafe , n S " ch a d»tectnc layer 
and Acoustical Microelectronics, pp. 223-247, Poly- =™ te * e comb,ned 8» to ™ de and fi f ld ° xide of a 
technic Press, New York, 1974. ^MOS dev.ce. The non-plananty of the dielectric layer may 
.^r, x,... , , . , , „ „ . ,, w j ■• alternatively be attnbuted to the non-plananty of underlying 

(11) F. H. Dill, J. A. Tut.le, A. R Neureuther, Modeling J5 , rench ^ mesfls of Qther topog ; apni / features 0 / lh f 

Positive Photoresist, Proceedings, Kodak Microelec- 

ironies Seminar, pp. 24-31, 1974. a polysilicon or other conductive layer is deposited 

(12) C. Mac, "Development of Positive Photoresists," conformably on top of the dielectric layer. An anti-reflective 
Journal of the Electrochemical Society, Vol 134, Janu- coating (ARC) may be optionally deposited on the conduc- 
ary 1987. 60 tive layer to reduce undesired reflections in a following 

(13) M. Stein, "Large Sample Properties of Simulations exposure step. 

using Latin Hypercube Sampling," Technometrics, Vol. The deposition of the conductive layer and optional ARC 

29, No. 2, pp. 143-151, May 1987. layer is followed by a spinning-on or other deposition of a 

(14) M. D. McKay and R. J. Beckman, "Using Variance photoresist (PR) layer. The PR layer may or may not be 
to Identify Important Inputs," 1994 ASA Proceedings, 65 planarized depending on process specifics. 

Section on Physical and Engineering Sciences, The photoresist-coated wafer is then positioned within a 

Toronto, August 1994. stepper by an alignment mechanism. Tiled areas of the 
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photoresist (PR) are successively exposed to a stepped this is accomplished by varying selected inputs or noise 

pattern of resist-modifying radiation. After the step-wise factors, then by running computer code to simulate the 

exposure, the photoresist is "developed" by, for example, variation likely to occur in manufacturing, 

soft baking the wafer to induce cross-polymerization in the Different process analysis and design approaches exist 

photoresist material and by subsequently dissolving away 5 with r tQ the descri tiorlj mo deling and modification 

non-hardened portions of the photoresist with a specific Qf a ^ M tfae variation {Q Qccm 

solvent. The developed photoresist defines an etch mask. . , & ^ c . /a. 

TTie developed wafer is next etched, for example in a m tne P rocess ' ™ e first a PP roach ' generally referred to as 

plasma etch chamber, so as to transfer the hardened image ^em design (as defined by Gemchi Taguchi in Introduc- 

in the photoresist to the conductive layer. This produces a 10 tion to Quality Control," p.75 (Asian Productivity 

corresponding pattern of lines in the conductive (e.g., doped Organization, 1986)), comprises the general product 

polysilicon) layer. architecture, schematic, functionality, and feature set. In the 

The photoresist mask is stripped off, or kept depending on second approach, generally referred to as parameter design, 

process specifics, and further process steps follow. One as defined by Taguchi, id. at p. 76, nominal process settings 

example of a further process step is the selective implant of 15 are determined in order to minimize performance variability 

dopants into exposed semiconductor regions so as to create and ^ In a mird approachj referred t0 ^ tolerance design 

self-aligned source and drain regions at opposite sides of ^ ^ 7g) ^ nmi[s are 

each conductive line, where the conductive line lies over u- ■ ■ u- a 

gate oxide. The width of the conductive line at such a region determined in order to achieve engineering objectives and 

of dopant implant defines the channel length of the formed 20 minimize costs. The invention described herein is in the 

IGFET transistor, tolerance design area. 

Within each of the above-described process steps, there Computer experiments have emerged as one way engi- 

are one or more variable physical attributes (or "process neers can model process designs using statistical methods to 

parameters") that control the final outcome of the produced devise process parameters for semiconductor manufactur- 

device. Some process parameters may be adjusted by a line 25 m g. The models allow engineers to increase product yields 

operator. Some are "set" by the design of the process by selecting variances in certain process factors such that the 

equipment that is installed into the mass-production line. V ariance(s) of the output characteristic(s) does(do) not 

Post-exposure development time and temperature are also ^ affect t fonnancCt 
subject to variance away from pre-established goal values. ^ . ,. j- ■ i a c *u 
The diffusion length of the development chemistry may vary 30 ° ne f od ^8 paradigm involves a description of the 
across a wafer. In a subsequent plasma etch, the variables cffect of a chan & e 10 one P rocess variatl0n on the out P ut of 
can include: time, pressure, temperature, flow rate, and another relative to some common reference. For the pur- 
field -proximity effects resulting from the pitch and step poses of this description, this paradigm will be defined as an 
profile of closely-spaced mask features. "engineering process window." In one aspect, a process 

Because successive steps of IC production tend to be 35 window is the inherent effect of a change in one parameter 

interdependent, a slight variation in parameters) of one resulting in an effect on any other parameter which may be 

process step can be magnified by a further variation in the plotted in two (or perhaps more) dimensions. Stated still 

parameters of a second process step to produce unacceptable anQther wayj ^ me ^ of twQ input parameterSj a process 

numbers of defective product at the output end of the ^ fa ^ description of thc tradeoff between the 

" to P e^k° tf^PR (photoresist) thicknesses decrease " tolerance of one input parameter and the tolerance of other 

slightly and the focal depth of the exposure optics also input parameters. Parameter design gives a process wmdow 

decreases slightly and the exposure dosage also decreases which defines, as a function of input parameters X r . . 

slightly during production of a first-sampled IC chip as the region in which the product performs in conformance to 

compared to the corresponding process parameters for a 45 specifications. Tolerance design yields, for a given range of 

second-sampled IC chip, the combined effect may be to input parameters, the tolerances X^Aj, . . . , X / ±A / which 

significantly shift the position and intensity of the radiation allow the product to perform in conformance to its perfor- 

exposure pattern relative to the photoresist layer during the mance specification. 

production of the first-sampled IC chip. The second-sampled i D computer experimentation involving engineering pro- 
IC chip may come off the production line in acceptable form 50 ceS s windows, calibration to realistic measurements is gen- 
while the first-sampled chip comes out of the same mass- era u y compared with computer-implemented models. Data 
production line in defective form. f rom the pr0C ess may be sampled from sacrificial areas of 

It is hard to pinpoint why mass-production yield for a actual processed substrates, with the sample data then used 
given circuit layout on a given mass-production line with empirical experiments or computer simulators (such as, 
becomes unacceptably low Tht statistical variance of PR 55 for & { ^ DEPIC T photolithography simulator avail- 
thickness across the production lot may be small. The aWc frQm Technol Modeling Associates of Palo Alto, 
statistical variance of focus across the production lot may be & & model Qf ^ & ^ 
small. But the physical interaction between the two noise n • j uu n * n * 
quantities can be such that the overall lithography process , P fOCess ^ows which compare results of altering two 
produces chips having a much larger variance in terms of 60 factors »" relativel y s * m P le: the model deals with changes 
critical dimensions. If the overall variance in final critical to one parameter which affect another parameter. The analy- 
dimensions becomes too large, production yield may suffer sis becomes much more complex where a number of input 
significantly. factors are involved as in, for example, a cross-correlation of 

The interdependence of variance of a number of specific device parameters which are utilized in a semiconductor 

process steps is not easily perceived. Hence, computer 65 fabrication process and affect conductive line formation, 

experimentation has aided the engineer by allowing distri- Such factors include, and are abbreviated in this specifica- 

butions of key output parameters to be estimated. Generally, tion as follows: 
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One descriptive method used, successfully to compute 
variance parameters is known as analysis of manufacturing 
variance (AMV) and is set forth in Heavlin, "Variance 
Components and Computer Experiments," reference paper 
no. 9, cited above. 

In analysis of manufacturing variance, the goal is to 
decompose the total variance of a process into components, 
one component associated with each input factor, and con- 
sistent with the variation to be anticipated in manufacturing. 
In a basic form of AMV, each noise factor is perfectly 
controlled, and the reduction in output variation noted. AMV 
allows for control of subsets of factors as well. 

AMV in its simplest form comprises using Latin hyper- 
cube (LHC) sampling for factors X v . . in a computer 
simulation, calculating the total variance of the resulting 
values, and substituting the mean of each factor to calculate 
the reduction in variance from the total variance to estimate 
the contributions of the factors held constant. AMV analysis 
is illustrated in the table shown in FIG. 3 A. In FIG. 3A, each 
X represents a noise factor for an input under determination 
where X«0, perfect control (no variance) of the input factor 
is presumed. One factor per iteration is presumed perfect, 
and all other factors are allowed to vary. 

Another methodology related to AMV is described in U.S. 
Pat. No 5,646,870 and referred to therein as incremental 
leveraging. Initially, a predefined goal for variation is 
selected — for example 10% — and an evaluation of each 
factor made on this basis. Process parameters are selected 
one at a time and the variance of the selected process 
parameter is temporarily reset to zero. The noisy set of 
simulations is re -run and the new total variance is calculated 
and stored while still retaining the variance value of an 
original first run where none of the relevant variances were 
reset to zero. 

Incremental leveraging has the advantage that it is good at 
suggesting one solution for input tolerance design. It does 
not, however, provide a mechanism for flexibility in allow- 
ing more than the one solution it provides for changing the 
tolerance constraints of the factors in the process under 
consideration. 

In co-pending application Sen No. 09/130,528, entitled 
Statistically Based Process Windows, inventor William D. 
Heavlin, assigned to the assignee of the present application, 
incremental leveraging methodology is utilized as a basis to 
define a statistically-based, process window algorithm 
wherein the response variability of a process is represented 
as a function of the tolerances of each variable. In brief, one 
utilizes the tolerance variations defined in the aforemen- 
tioned algorithms as inputs to determine, for any number of 



process variables, the effect of varying the tolerance of each 
factor on the tolerances of other factors under consideration 
relative to a response factor for the entire process. This 
algorithm may be performed in a computer and have an 

5 output provided to a process engineer who may then deter- 
mine optimal modifications to the nominal process under 
consideration to derive a toleranced process. 

In general, the method defined in co-pending application 
Ser. No. 09/130,528 comprises: defining a nominal process; 

io evaluating F input factors {X^f=l, . . . , F} which affect the 
output of the nominal process; and defining, as a result of the 
evaluation, a toleranced process having an output response 
less susceptible to variance in the input factors. The evalu- 
ation is performed by representing the variability of the 

is response of the process as a function of the factors' current 
variability and selected tightening factors; performing an 
evaluation where the tightening factors are varied using an 
orthogonal array, varying the factors over [0,1]^; and 
modeling, using kriging interpolation, a neural network or 

20 equivalents, to determine moments or other statistics of a 
distribution of the process for given changes in the con- 
straints of each factor. 

The method of forming the semiconductor device in 
accordance with the present invention results in better 

25 planning of next-generation equipment requirements and 
ultimately in improved product yields as the considered 
factors are improved upon. 

SUMMARY OF THE INVENTION 

30 An objective of the invention is to allow process engi- 
neers to determine, on a lot-by-lot basis, whether and how to 
make improvements in the processing of the particular lot 
under consideration. 

A further objective of the present invention is to provide 

35 an evaluation method for process engineers to make quali- 
tative judgements about the manufacturability of certain 
processes based on models of the output given a variety of 
inputs to the process. 

Yet another object of the invention is to improve yields in 

40 multi-faceted manufacturing processes, and particularly 
semiconductor manufacturing applications. 

In one aspect, the invention, roughly described, comprises 
a method of fabricating a multi-component semiconductor 
device assembly, comprising: defining a nominal semicon- 

45 ductor manufacturing process having a plurality of process 
steps, the nominal process including a number of process 
input factors X^Xp which, when perfectly controlled to 
meet the nominal process settings for such factors, produce 
a device assembly with defined operating characteristics: 

50 performing at least one of said plurality of process steps on 
at least one lot of wafers; and evaluating the process by: (1) 
representing the variability of the response of the operating 
characteristic as a function of the actual variability of at least 
one of said factors and at least one tightening factor; (2) 

55 evaluating said input factors subject to a variety of tighten- 
ing factors array over [0,1 J^; and (3) modeling the output 
using an interpolation function to determine suitable toler- 
ance models impacting control of one or more of said input 
factors. 

60 BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be described with respect to the 
particular embodiments thereof. Other objects, features, and 
advantages of the invention will become apparent with 
65 reference to the specification and drawings in which: 

FIG. 1 is a flow diagram illustrating one embodiment of 
the method of the present invention. 



05/28/2004, EAST version: 1.4.1 



DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 



US 6,708,073 Bl 

7 8 

FIG. 2 is a flowchart illustrating the evaluation portion of mean and a deviation about the mean, and also between 

the present invention. manufacturing-capability -defined lower and upper limits, 

FIGS. 3A and 3B illustrate the input and output, e.g., mean ±3 standard deviations, 

respectively, of an incremental averaging algorithm process In accordance with the invention, at some point in the 

over 210 intervals. 5 series of steps 10-14 performed on a lot of wafers, an 

FIG. 4 is an interactive dialog table illustrating a portion analysis of the lot-specific factors is performed and a deci- 

of the method of the present invention. sion regarding the wafers in the lot can be made resulting in 

FIG. 5 is a flow chart illustrating the statistical process either J ^vings with respect to further processing or a 

window evaluation portion of the method of the present in ^working of the ^lot prior to continued processing, to correct 

+ . iu tne problematic factor, and improve the resulting yield ot the 

invention. . C, . ' . r . . , , & J . , 

lot. This contrasts with previous statistical analysis wherem 

HG. 6 is a graphical representation of the command the entire process the bt was tQ ^ subjected t0 might 5e 

instructions and matrix input for the method set forth in FIG. analyzed prk)r tQ processing< M shown m Fia lf each 

5 " process step 10-14 results in a given set of data 10a-14a 

is which has a given standard deviation 10b-14b for the lot. As 
such, the practical application of the statistical analysis of 
the performance factors allows for decisions to be made in 

In manufacturing processes, such as those utilized for real time in the processing of a given lot of wafers, 

manufacturing a semiconductor device, imperfect control of In FIG. 1, three processing steps are completed on a given 

input variables propagates through the manufacturing pro- 20 lot when analysis in accordance with the present invention 

cess to induce variation in the manufactured device, and in occurs. However, it should be recognized that the analysis of 

response variables describing the device. Where there are a given data may be performed after any number of process 

number of process inputs, one can define the relative impor- steps 1-N where N is the last step in the process. This further 

tance of such inputs, and more than one set of input differs from the statistical analysis of co-pending application 

tolerances which will meet the desired response. This means 25 Ser. No. 09/130,528 in that it substitutes lot-specific inputs 

that there are trade-offs in the control tolerances, or required for computed inputs to adapt the analysis to the particular 

distributions of standard deviations, among the various input lot. 

factors. While process windows represent one-to-one trade- The evaluation step 100 of the present invention may be 

offs, a more generalized model is required to consider all performed in an independent computer unit or can be a 

sets of input tolerances. The instant application presents a 30 module within a computer that also carries out functions of 

novel tolerance design method for manufacturing semicon- other described or to-be-described modules. The output of 

ductor devices. the model may be used by the process engineer in possession 

While the method of the present invention will be of other data, such as cost factors, to define constraints and 

described herein with respect to a semiconductor manufac- 35 make choices for altering the semiconductor manufacturing 

luring process, it should be recognized that the method of the process to determine which constraints to confine or relax, 

present invention may be applied to any process wherein a Following the evaluation step 100, the decisions to scrap 

number of process variables affect the output variance of the 22, rework the lot 24, send the lot to higher capacity 

manufactured product. equipment 26 or take some other action 28, can be made 

FIG. 1 illustrates the method of the present invention. 40 based on lot-specific data to result in an overall improved 

FIG. 2 is a block diagram representing the basic statistical yield for the lot. 

process window methodology as utilized in the present FIG. 2 is a flowchart depicting, in one aspect, the evalu- 

application. ation method in accordance with the present invention. At 

FIG. 1 shows a series of process steps 10-14 which are step 55, a manufacturing process under consideration is 

performed to manufacture a typical semiconductor device. 45 defined based on the particular characteristics of the product 

By way of example, a typical semiconductor mass- to be manufactured. In the case of a semiconductor device, 

production line which is part of an overall wafer processing this process will include a number of process steps as 

system may include: (1) a station for coating each wafer with described herein. The process defined to produce the product 

material layers such as the anti-reflective coating and pho- is generally referred to as the "nominal" process and results 

toresist layers; (2) a station for exposing the material layers 50 in process steps, each process step having one or more input 

(the photoresist layer) to a radiation pattern; (3) a station for factors, such as those factors set forth above with respect to 

developing each exposed wafer; and (4) a station for etching the semiconductor manufacturing process, 

each developed wafer. The method of the present invention As shown at step 60, for a set of process steps, the input 

allows tweaking of the process at any one or multiples of factors whose variance affects the output yield of the 

such stations. Steps 10-14 represent individual steps at any 55 process, based on their interaction with other factors, are 

one of such stations. defined. Each input factor will have a nominal or mean 

Typically, each post-etch wafer lot is moved further value, and a standard deviation about the mean value. The 

downstream for further processing (e.g., dopant implant) in standard deviation can be an actual deviation during the 

accordance with well-known techniques. Before such fur- process or a simulated deviation for purposes of the evalu- 

ther processing, sample wafers from each post-etch wafer lot 60 ation of the present invention. 

may be sampled to determine critical dimensions in order to Step 100, defined in dashed lines and comprising several 
calibrate the computer model of the method of the present sub-steps, comprises the evaluation portion of the method of 
invention. As shown in FIG. 1, step 100 represents the the present invention. The evaluation portion 100 can corn- 
statistical process window evaluation used over a set of F prise a set of instructions which are processed by a 
individual factors over N steps. For each component, the 65 computer, in one or more modules as described above, 
simulation is used to generate a statistical distribution of designed specifically to implement the evaluation based on 
each fit-determining parameter. The distributions have a the input factors provided to it. 
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In one embodiment, as depicted in FIG. 2, the evaluation 
portion 100 of the method of the present invention 
comprises, at step 110, defining the variance of the manu- 
facturing output — the product produced — as a function of 
the variability of the input factors and tightening factors 5 
applied to each input factor. Once so defined, at step 120, an 
evaluation is performed by varying the tightening factors to 
determine alternative input tolerances which yield accept- 
able manufacturing outputs. Finally, at step 130, the evalu- 
ation is modeled to determine the process variance as a JQ 
function of the tightening factors. Following this evaluation, 
a determination as to the fate of the lot can be made with an 
intelligent projection of the performance of the lot. 

FIGS. 3 A and 3B illustrate the output of the incremental 
leveraging calculation for a given simulation. This output of 
the incremental leveraging algorithm provides the tightening 15 
factor starting points for determining the matrix of tighten- 
ing factors used in the evaluation portion of the present 
invention. The input table 40 includes the nominal mean and 
the standard deviation values for seven input factors con- 
sidered in this example which comprises line width defini- 20 
tion process, similar to that described above and in U.S. Pat. 
No. 5,646,870. Each factor is assigned an effort of 0.9, 
reflecting that, for this exercise, the design goal is a 10% 
incremental change. It should be recognized that other effort 
factors may be utilized in accordance with the present 25 
invention. The factor of a 10% increment is selected based 
upon common practice in the industry. The output table 50 
shown in FIG. 6B comprises a calculation in accordance 
with the incremental leveraging algorithm described with 
respect to FIGS. 5A-5B, over 210 increments. Column 1 of 30 
output table 50 shows the resulting weighting of weighted 
tolerance factor for each of the seven factors analyzed. For 
example, for incremental row 205, the focus factor has a 
standard deviation of 0.06960, and the factors of exposure 
(ExpD), focus and post exposure bake (PEB) are shown to 35 
have the greatest contribution over these increments. Incre- 
mental leveraging assigns the focus parameter 7 increments, 
exposure 11 increments, and PEB 12 increments. This yields 
tightening factors for ExpD of 0.31381=(0.9) u , for Focus of 
0.47830=(0.9) 7 , and for PEB of 0.28243<0.9) 12 . 40 

FIG. 4 shows an interactive dialog for a statistical process 
window algorithm that is used in the evaluation portion of 
the present invention. FIG. 5 presents a generalized descrip- 
tion of the evaluation algorithm of the present invention 
which should aid in understanding the more specific descrip- 45 
tion to follow. In general, the dialog depicts the interaction 
between terms necessary to loosen the tightening factor for 
focus from an interactive leve raging-derived value to 1.0, 
and the effect on the tightening factors of all other factors 
under evaluation. It should be noted that while focus is 50 
illustrated as the factor under consideration, other factors, or 
multiple factors, may be considered. 

The tightening factors for exposure, focus and post 
expose bake are calculated tightening factors based on an 
incremental leveraging algorithm, as described above. All 55 
other factors are assigned a factor of 1.0=(0.9) 0 . A resulting 
control factor is computed. The control factor reflects the 
3 a/mean value as computed in an incremental leveraging 
algorithm. Next, the tightening factor for focus is relaxed 
(from 0.47830 to 0.081 (at col. 2)), and all other factors 60 
except ExpD and PEB tightened, and the resulting control 
recalculated. As will be noted, the control factor increases to 
0.11536. Next all factors, including ExpD and PEB are 
tightened, and the control factor again calculated. This 
process continues as illustrated in FIG. 4. 65 

The goal of the aforementioned dialog can be seen by 
reviewing the subsequent rows 4-11 wherein the focus 
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contribution is relaxed to zero increments (0.9)°«1.0 in row 
6, as all other factors are tightened until the resulting control 
factor achieves proximity with the initial, incrementally- 
dcrived control factor in row 1 of the dialog. 

Comparing the statistical process window dialog with 
incremental leveraging, incremental leveraging gives an 
output for each factor of (a)^)*^, a non-negative integer, 
for each tolerance. With SPW, the resulting control allows 
one to relax or modify any number of the factors (not merely 
the single focus factor identified in the dialog of FIG. 4), to 
achieve the resulting control value. It should be noted that in 
the method of the present invention, each of the process 
factors may be varied relative to external considerations, 
such as cost or ease of control over the factor. In the 
aforementioned example of FIG. 4, the noise factor is 
selected at (0.9). In one aspect of the method of the present 
invention, the noise factor is varied over [0,1] F 
(mathematically, the interval from 0 to 1, inclusive) in an 
orthogonal array by means of a computer simulation as 
described below. In an alternative aspect, the noise factor 
may be varied over [0 t CY, where C^l, to incorporate both 
tightening and loosening factors into the evaluation. 

The evaluation portion 100 of the method of the present 
invention will be described with reference to FIGS. 5 and 6, 
FIG. 5 is a flow chart of the evaluation algorithm 100 of the 
present invention. FIG. 6 is a representation of the command 
instructions and mathematical inputs of the evaluation 
method of the present invention. 

FIGS. 5 and 6 consider, for example, a lithography portion 
of a semiconductor manufacturing process. A nominal 
lithography process has been determined by process 
engineers, based on the device to be manufactured, and the 
nominal process will be evaluated in accordance with the 
inventive method. 

As shown in FIG. 5, process step 110a consists of defining 
a resulting linewidth critical dimension as a mathematical 
function of a number (f) of input factors. The result is an 
output variable (cd_dc5a) which is the goal of the linewidth 
definition process (e.g., a linewidth). In one aspect of the 
present invention, this definition step 110a can be performed 
by a kriging interpolation in accordance with the approach 
outlined in B. D. Ripley, "Spatial Statistics," article 7 set 
forth above, or by suitable fast computer simulation tech- 
niques. Essentially, step 110a can be decomposed into the 
substeps of: building; a kriging model to determine the 
response dimension relative to the input factors; and cali- 
brating the model. The calibration adjustment may be per- 
formed in accordance with the teachings of W. D. Heavlin 
and L. Capodieci in "Calibration and Computer 
Experiments," cited as article 15, above, wherein simulated 
critical dimensions are compared against empirical critical 
dimensions. The calibration model may result from either 
kriging interpolation or use of a neural network in accor- 
dance with the teachings of the aforementioned reference. 
The output critical dimension used in the method of the 
present invention (cd_dc5a) is the result of the calibration 
model and the simulation critical dimension resulting from 
the kriging model. 

In FIG. 6, the evaluation steps 110a, 115 and 120 a are 
represented as a single command line 90 which may be used 
to call the functions or tables described herein. The tables 
and arrays described with respect to FIG. 5 are graphically 
represented. In FIG. 6, at step 115, the nominal process 
values of the input factors and the starting process tolerances 
are defined for input to the evaluation portion of the method. 
The input factors (mean) and variances (stdev) are set forth 
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in a table (pw__mst). The nominal process targets at the input 
factors are the starting (wide-standard deviation, current 
generation) process tolerances of the input factors. 

In accordance with the present invention, the input factors 
for the mean and variances for each factor may be real data 5 
10a-14a, or calculated data in accordance with the forego- 
ing description. 

The current generation values of the input factors are then 
used to build a data set for the evaluation-portion of the 
method of the present invention. 

At step 120a, a schedule of tightening factors (L49 as 
illustrated in FIG. 6) is then evaluated against the defined 
output variable (cd_dc5a), which is subjected to a Latin 
hypercube sampling (LHS343) for each row of the schedule 
of tightening factors in accordance with the following 
pseudocode: 



loop "schedule" (L49) 
for row i of schedule 
copy LHSj from LHS 
loopj 

Col j of LHS i -mst(mean)(j)+(Col j of LHSJ* 
mst(std)(j)*sched(ij) 
end loop j 
LHS; (343 rows) — cd_dc5a 
m(i) - mean of LHS 343 -* cd_dc5a 
s(i) - stdev of LHS 343 -*■ cd_dc5a 
end loop i 



25 



30 



Step 128 results in an output (pw__dc5cd as represented in 
FIG. 6), which is a mean and standard deviation table with 
a number of rows (49) corresponding to the number of rows 
(49) in the schedule (L49) of tightening factors. The output 
pw_dc5cd is then used in a kriging interpolation at step 35 
130a, to yield a mathematical function which is the standard 
deviation of the output response and the tightening factors. 

The output of the kriging interpolation shall be referred to 
herein as the "statistical process windows (SPW) function." 
The SPW function is the standard deviation of the output 40 
response as a function of each input factor and a correspond- 
ing tightening factor for the particular input factor. 

As represented in FIG. 6, at step 94, the log of col. 1 
(stdev) of pw_dc5cd is taken to exaggerate small values of 
the output variance. As shown in FIG. 6, in the particular 45 
kriging estimation which is utilized in the present 
invention — kpw_dc5cd, two iterations are used with a ridge 
(a.k.a. nugget) parameter of 0.01. 

At step 98 of FIG. 6, the exponent is taken in accordance 
with conventional statistical techniques to invert the log 50 
function of step 94. 

As shown at step 98, the kriging interpolated estimate is 
a deviation of 0.060705, compared to a directly-calculated 
value of 0.060533 for one particular point of interest. 5S 

Once the evaluation portion of the present invention is 
completed for a particular step, the process evaluator can 
make any number of decisions as indicated in FIG. 1. 

The many features and advantages of the present inven- 
tion will be readily apparent to one of average skill in the art. 60 
The method of the present invention provides a model which 
can be used in processing devices based upon computer- 
simulated results of processes, or on empirically derived 
models of each process. No specification of a cost-of- 
tolerance function is required, thereby simplifying applica- 65 
tion of the method to various processes or sub-processes in 
device manufacturing. 



The method of the present invention utilizes an evaluation 
process which makes explicit the trade-off in the control 
requirements amongst various manufacturing steps, and 
provides a quantitative analysis of this trade-off for use by 
process engineers in developing the fabrication process. 
Process engineers are provided with an explicit role for 
engineering judgment by the quantitative description of each 
trade-off in the process. 

Obvious modifications to the method of the invention will 
be apparent to those of average skill in the art. For example 
and without limitation, at step 130a, other methods of 
modeling the output of step 120a, such as utilizing a neural 
network, are contemplated. All such modifications are 
intended to be within the scope of the invention as defined 
by the claims. 

What is claimed is: 

1. A method of fabricating a multi-component semicon- 
ductor device assembly, comprising: 

(A) defining a nominal semiconductor manufacturing 
process having a plurality of process steps, the nominal 
process including a number of process input factors 
Xj-Xf which, when perfectly controlled to meet the 
nominal process settings for such factors, produce a 
device assembly with defined operating characteristics, 
the input factors each having a distribution of accept- 
able tolerances which are acceptable inputs to the 
process which result in the defined operating charac- 
teristics; 

(B) performing at least one of said plurality of process 
steps on at least one lot of wafers; 

(C) evaluating the process by: 

(1) representing the distribution of the response of the 
operating characteristic as a function of the actual 
distribution of at least one of said factors and param- 
eterized by at least one parameter; 

(2) evaluating said input factors subject to a variety of 
tightening or loosening factors in an array over 
[0,Cf , where C^l; 

(3) modeling an output using an interpolation function 
to determine suitable tolerance models impacting 
control of one or more of said input factors. 

2. The method of claim 1 wherein said step C(3) com- 
prises evaluating said factors in an array over [0,1 J*. 

3. The method of claim 1 wherein said step C(3) com- 
prises evaluating said input factors in an array over [0, Cf 
where Ci^l. 

4. The method of claim 1 wherein said step (B)(1) 
comprises the sub-steps of: 

defining an output response based on a matrix of said 
factors; and 

calibrating said output response based on empirical data 
of the output response. 

5. The method of claim 4 wherein said sub-step of 
defining comprises utilizing a kriging interpolation. 

6. The method of claim 4 wherein said sub-step of 
defining comprises utilizing a neural network. 

7. The method of claim 1 wherein said step (B)(2) 
comprises: 

for a schedule of tightening factors, evaluating for each 
row of the schedule, a Latin hypercube sample of, for 
each column of the sample, at least the mean and 
standard deviation of each input factor, taking the mean 
of the Latin hypercube sample output of each row of the 
schedule in the standard deviation of the Latin hyper- 
cube sample for each row of the schedule. 

8. The method of claim 1 wherein said step (B)(3) 
comprises applying a kriging interpolation to the output of 
said step (BX2). 
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9. The method of claim 1 wherein said step (B)(3) 20. The method of claim 13 wherein steps (C)(2) and 
comprises utilizing a neural network on the output of step (C)(3) are repeated based on empirical choices made accord- 

(B) (2). ing to external factors. 

10. The method of claim 1 wherein steps (B)(2) and (B)(3) 21. The method of claim 13 wherein said step (C)(1) 
arc repeated based on empirical choices made according to 5 comprises: 

external factors of the nominal process. simulating an output response as a function of the input 

11. The method of claim 1 wherein said step (B)(1) factors to determine an output variable; and 
comprises: calibrating the simulation by comparing empirically 

simulating an output response as a function of the input determined output responses as a function of selected 

factors to determine an output variable; and 3Q factors to simulated output responses generated by the 

calibrating the simulation by comparing empirically simulation, 

determined output responses as a function of selected 22. The method of claim 21 wherein the step of simulating 

factors to simulated output responses generated by the comprises using a kriging interpolation, 

simulation. 23. A method of evaluating a lot of wafers being pro- 

12. The method of claim 11 wherein the step of simulating cessed during a series of steps in a nominal manufacturing 
comprises using a kriging interpolation. process having a plurality of components, each component 

13. A method of fabricating a multi-component assembly, having a mean and standard deviation; comprising the steps 
the method including at least components Xl-Xf wherein of: 

each of said components Xl-Xf has a distribution of accept- (A) performing at least one of said series of steps in said 

able tolerances within which the method will produce the manufacturing process to determine an actual current 

multi-component assembly in an acceptable fashion, com- 20 variability for a factor; 

prising: (B) representing the variability of a response of a system 

(A) defining nominal values of each component; to the current variability of at least one of said factors 

(B) performing at least one of a series of steps on a lot of and at lfe ast one tightening factor; 

wafers to determine a real variance of at least one of ( c ) evaluating said components against a schedule of 

said components; 25 tightening factors in an orthogonal array; 

(C) evaluating the performance of the method on the lot (D) modeling an output of the evaluation using interpo- 
5y. lation to determine a mathematical function defining a 

(1) defining the distribution of at least one critical toleranced process; and 

dimension of the assembly as a function -of said 3Q (E) applying at least one of the models to a fabrication 

actual distribution of at least one of the components system to complete the model, 

and parameterized by at least one parameter, at least 24. The method of claim 23 wherein said step (A) 

one parameter associated with each component; comprises: 

(2) analyzing the nominal components by comparing a defining an output response based on a matrix of said 
Latin hypercube sampling of the nominal compo- 35 factors; and 

nents against a schedule of parameters; calibrating said output response based on empirical data 

(3) modeling the analyzed components to define the of the output response. 

variance of an output as a function of the compo- 25. The method of claim 24 wherein said step of defining 

nents and the parameters; and comprises utilizing a kriging interpolation. 

(D) manufacturing the multi-component assembly by 40 26. The method of claim 24 wherein said step of defining 
performing the method and controlling ones of said comprises utilizing a neural network. 

components based on the output of step (B). 27. The method of claim 23 wherein said step (B) 

14. The method of claim 13 wherein said step (C)(1) comprises: 

comprises the sub-steps of: for a schedule of tightening factors, evaluating for each 

defining an output response based on a matrix of said 45 row of the schedule, a Latin hypercube sample of, for 

factors; and each column of the sample, the mean and standard 

calibrating said output response based on empirical data deviation of each input factor, taking the mean of the 

of the output response. Latin hypercube sample output of each row of the 

15. The method of claim 14 wherein said sub-step of schedule in the standard deviation of the Latin hyper- 
defining comprises utilizing a kriging interpolation. 50 cube sample for each row of the schedule. 

16. The method of claim 14 wherein said sub-step of 28. The method of claim 23 wherein said step (C) 
defining comprises utilizing a neural network. comprises applying a kriging interpolation to the output of 

17. The method of claim 13 wherein said step (C)(2) said step (B). 

comprises: 29. The method of claim 23 wherein said step (C) 

for a schedule of parameters, evaluating for each row of 55 comprises utilizing a neural network on the output of step 

the schedule, a Latin hypercube sample of, for each (B). 

column of the sample, the mean and standard deviation 30. The method of claim 23 wherein steps (B) and (C) are 

of each input factor, taking the mean of the Latin repeated based on empirical choices made by external 

hypercube sample output of each row of the schedule in factors. 

the standard deviation of the Latin hypercube sample 60 31. The method of claim 23 wherein said step (A) 

for each row of the schedule, comprises: 

18. The method of claim 13 wherein said step (C)(3) simulating an output response as a function of the input 
comprises applying a kriging interpolation to the output of factors to determine an output variable; and 

said step (C)(2). calibrating the simulation by comparing empirically 

19. The method of claim 13 wherein said step (C)(3) 65 determined output responses as a function of selected 
comprises utilizing a neural network on the output of step factors to simulated output responses generated by the 

(C) (2). simulation. 
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32. A computer program embodied in a non-volatile 
memory, responsive to inputs, comprising: 

(A) means for defining an output factor based on the 
variance of the inputs; 

(B) means for evaluating input factors by computing: 
for each row i of a schedule of input factors, 

for each column in a Latin hypercube sampling 

Col j of L7/5 r mst(mean)0>(Col / of I/^S / )*mst(std)0')*sched(t,/) 
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end loop j, to compute 
m(i)«mean of LHS rows 
s(i)=stdev of LHS rows; and 

5 (C) modeling the tolerance of the input factors by kriging 
interpolation to yield an output variance which is a 
function of the inputs. 

* * * * * 
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